Essential Role of the Cooperative Lattice Distortion in the Charge, Orbital and Spin 

Ordering in doped Manganites 
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The role of lattice distortion in the charge, orbital and spin ordering in half doped manganites 
has been investigated. For fixed magnetic ordering, we show that the cooperative lattice distortion 
stabilize the experimentally observed ordering even when the strong on-site electronic correlation 
is taken into account. Furthermore, without invoking the magnetic interactions, the cooperative 
lattice distortion alone may lead to the correct charge and orbital ordering including the charge 
stacking effect, and the magnetic ordering can be the consequence of such a charge and orbital 
ordering. We propose that the cooperative nature of the lattice distortion is essential to understand 
the complicated charge, orbital and spin ordering observed in doped manganites. 

PACS numbers: 75.30.Vn, 7f .45.Lr, 75. 25. -Hz, 7f .70.Ej 
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Pri/2C^i/2Mn03lll, Lai/2Sr3/2Mn04BB and 



Lai/3Ca2/3MnO 



The unusual charge, orbital and spin order- 
mg (COSO) irL. manganites have recently attracted 

below a certain temperatures Tco, 
electronic carriers become localized onto specific sites, 
which display long-range order throughout the crystal 
structure (charge ordering). Meanwhile the filled Mn'^^ 
Cg orbitals also develop long-range order (orbital or- 
dering). When the temperature is further decreased 
to a much lower temperature Tjv, an antiferromagnetic 
(AF) magnetic ordering with a zigzag pattern sets in 
(Fig. |l|(a)). In some others like Lai/2Cai/2Mn03El and 
Ndi/2Sri/2Mn03l!l, the system first undergoes a ferro- 
magnetic (FM) transition at the Curie temperature Tc, 
then enters into the CE-type charge, orbital and spin or- 
dering (COSO) state at a lower temperature Too — Tn- 
Theoretically, it has been proposed that the charge and 
orbital ordering (CQQ.) in half-doped manganites has a 
magnetic spin originoEI. However, such a theory can not 
be applied to those materials with Tco > Tn, where the 
COO is established before the spin ordering. Even for 
those whose Tco = Tn, it has been shown that in a pure 
electronic model, the on-site repulsion U destabilizes the 
CE structure towards the rod-type (C-type\_AF state 
1(b)) in realistic parameter regime of [/El. Yunoki 
considered the effects of lattice distortions (LD), 
and found that both the noncooperative LD (NLD) and 
cooperative LD (CLD) can lead to the CE-type COSO. 
In their work U is neglected and the calculation was per- 
formed on a 4 X 4 X 2 lattice where the size effect is 
prominent. Since in the. absence of U the CE state can 
be obtained without LDB'lI, the role of LD seems not very 
clear there. It is desirable to clarify what the LD results 
for the CE state would be destabilized by U. 

In this work, we investigate the role of LD and large U 
in the COSO in half doped manganites. For fixed mag- 
netic ordering, by studying the competition between the 
CE and C states, we found that only CLD can stabi- 
lize the CE state in the presence of large U. Further- 



more, without invoking the magnetic interactions, the 
CLD alone can lead to the experimentally observed COO, 
including the charge stacking effect. The magnetic order- 
ing is the consequence of such a COO. 
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FIG. 1. View of the (a) CE and (b) C phase in the x-y 
plane. The arrows refers to the spin. Along the z-direction 
neighboring sites have the same charge and orbital states but 
opposite spins. 

The interaction concerning the lattice freedom includes 
two parts: Hiat — H^p -\- Heia- H^p is the coupli.ngj be- 



tween lattice distortion and eg electrons, given bylljil 
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where Ux, Oz are Pauli matrices, (£[ is electron operator of 
orbital 7 (2 = d322„^2,z = dx'^~y^\ ^.nd Qu — -^(wf + 
vl + 0, O2. = ^^{vl - vt),Q^., = ^^{2vl - V- - v\) are 
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the breathing (Qu) and Jahn- Teller {Q2i, Qsi) rnodes of 
the LD. Here v" = — with being the a- 

component of the displacements from the equilibrium po- 
sition of the neighboring oxygen ion in the^d^a direction. 
The parameter (3 is expected close to l/2Ea. The index 
of spin has been omitted, which is always parallel to the 
local spin due to the strong Hund's coupling. Through- 
out the paper, we use a = x,y or z to denote either the 
direction or the orbital state, in the latter case it refers 
to the orbital d^a^-r^, whose orthogonal state is denoted 
as a. There are relationships q'^ = cf/2 =F -\/3cf/2 and 
c^'^ = ±-\/3cf /2-|-cf/2. Heia is the elastic energy and de- 
pends on the relative displacements of neighboring atoms 
with respect to the ideal perovskitc lattice. In a unit cell 
of the perovskite Ai_2,A^Mn03, there are three kinds of 
atoms: Mn, O and Z(=A or A'). The main contribution 
to Heia may include the elastic energies of the neighbor- 
ing Mn-0, 0-Z and Mn-Z atoms. Up to now no studies 
concerning the elastic energy of Z atoms has been made, 
however, this energy should be important, as without it 
the Z atoms can have arbitrary displacement instead of 
sitting in the center of the cubic cell cornered by eight Mn 
ions, and experimental observatians indicate that the Z 
atoms also participate in the LDO. In harmonic approx- 
imation, the elastic energy of the CLD may be written 
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The tensor G = 
where P,D,R and W are 3x3 matri- 
= S^^,{3Ki + AK^)Sl/{-iKiSl + 4if3), 

C,,Wxx = 3KlhxC,Wxy = 

4i^3(l - C^Hx - ClSlyHy - 

y'^xz-^-'-zi I 'iK2{2 — Cxy — Sxy — Cxz ~ Sxz)iRxy = 

CxySxy[AK:i{ClHx + ClHy - SlH,) + QK2], where 5„ = 

si\l{qa/2) and Ca = COs((7q/2), Saa' = SaSa' and Caa' = 

CaC„s C = CxCyCz, = AK^/iSKiSl + AK3) and 
ha = HaSaCa- The othcr elements oi D,R and W can 
be obtained by exchanging the indices, e.g., Dyy 



and Dyx = 

nected to G through u; 



SxyhyC z 



hyC 



and 



etc. The displacements are con- 



fAtan(g„/2)(/">/i^i 



with Uq ct be- 
ing the Fourier transform of Ui+a, and Aq is a func- 
tion of Uq.Q, and ^q. It should be pointed out that 
the form of Eq.(||) is actually general for Hiat with 
any harmonic H^ia^ and different choice of H^ia leads 
to different G. 



For example, the H^ia in Ref. 13 



eludes Ki and K'l (the spring constant between neigh- 
boring Mn sites) terms, where the G tensor is ci^^, (q)_= 



5aa'{Ki + 2K[Sl)l{Ki + 2K[). While in Refs . |10| and |lj 



and G(f),(q) 



NLD and CLD yield G^^, = S^a 
respectively. Here the difference between CLD and NLD 
is whether the G tensor depends on q or not. 
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where 5, u and A are the dimensionless displacements 
of the Mn, O and Z ions with reference to the ideal 
perovskite lattice, e^, and e,, arc unit vectors along 
the directions of neighboring Mn-0, Z-0 and Z-Mn, re- 
spectively, with K,^ and 77 the indices of neighbors. In 
principle, the spring constants between Z-0 and Z-Mn 
depend on whether Z=A or A', here to simplify our 
study we replace them by the averaged K2 and K^. 
Since the distance between these neighboring atoms are 
Lyin-O < io-z < Lz-M-a, One expects that the spring 
constants Ki> K2> K^- nn 

In the classical treatment of the LDE2rE3, the displace- 
ments of various sites are determined by minimizing the 
total energy of the system, dHiat/ dwf = 0, where wf 
(w = (5, A or u) is the a-component of the displacement. 
For harmonic H^ia, such a calculation can be easily per- 
formed in the momentum space to get the optimized val- 
ues of the displacements. Then after substituting these 
displacements into Eqs.(|^) and (||), Hiat reduces to 
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where e/ = X^/Ki, /" is the Fourier transform of 
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FIG. 2. Calculated go, 91,92 and g'l as a function of 
K2/K1, with fixed K3/K2 = 0.5. 

Eq.@ indicates that LD results in an effective elec- 
tronic interaction. In real space, it is 
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eff ^ 
lat 



\ Z^WoJi h ~ 91 J I Ji+a - 92 J I Ji+2& 



+9[ J2 /r'/r;j+ E'^l'/"//}' (4) 

a'(^ci) ij,aa' 

where the sum includes all the other terms. g{ is 
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found to be smaller than the coupling coefficients of the 
first several ftfi+na terms, but larger than any other 
coefficients G]^^,. Fig. || shows the calculated values 
of (70,51,52 and g[. From Eq.(^, the main effects of 
the CLD corresponds to an effective short-range orbital- 
dependent coupling between occupied Mn sites. If G is 
replaced by G^'^ (i=l,2,3), then in the CLD cases of G^^^ 
and G^'^\ there are only the go and gi terms, while in the 
NLD case of G^^\ there is only the go term. 

Now let us investigate the effect of LD to the COSO 
at half doping. First we see the case with fixed CE-type 
spin ordering. In the one-dimensional zigzag FM chain 
(see Fig. 0(a)), the double-exchange (DE) Hamiltonian 
reduces to 

Hde = ^ [d^j(ticc,i+i + t2Cc,i~^i) + H.c] 
i=even 

+U{nB^n'B.^ + n^ca+inc,t+i)^ (5) 

where dsi — cf and cf for i — 4j and 4j -I- 2 (j is an 
integer), ca = (cf,±c|)-^ for i=4j -I- 1 and 4j -|- 3, and 
i =even and odd corresponds to the bridge (Mn'^+) and 
corner (Mn'*+) sites, respectively, ti^2 — ~{t/'2,^V^t/2), 
nsi = d^Bi'^Bi and n'^j = dg^d'^i, with dg^ being the 
orthogonal state of orbital dsi, and n^^j = Cq^Cq^. In 
realistic manganites, the parameter regime of the on-site 
repulsion [/ = 10 ~ 20t, so that the system is strong 
correlated. In this regime double occupancy of electrons 
at a site is almost forbidden, and it is appropriate to 
use the|-|Gutzwiller projection (GP) method, valid for 
U — > ooci, to take into account such a strong correlation 
effect. In GP we introduce constrained electrons at each 
site. Each electron operator in Eq.(||) is replaced by the 
corresponding projected operator to eliminate the double 
occupancy, — > (1 - ng^)dBi and c^^ ^ (1 - n])c1.^ 
(7 = z(z) when 7 — z{z)). Then we use a mean-field ap- 
proximation by decoupling the high order terms such as 

ter a Fourier transform, we have 

+ J2''4kd'Bk+E°DE, (6) 

fc 

where tj, = (1 — ^ "'c)^fc with tf. = —tcosk 

and tl = iVSisinfc, = -2(1 - n'B)ReJ2k{d^Bk^k'^ck)' 
e' = -2ReEfe^(l - 4)(4fcifc4fc). ^nd E^^ is the MF 
energy constant. 

For Hiat of Eq. (||) , in the case of CE type spin ordering 
the sum over q includes q = (0, 0, 0), ±{tt/2, 7r/2, 0) and 
(tt, tt, 0), and can be denoted as g = 0, ±7r/2 and tt in the 
ID FM chain. By decoupling the quartic term fg^fg 
{f^^)ff +f"Hff). Eq.(|) reduces to 

HL7 = + H.c.) + El„ (7) 

k,q 



where Ck is the Fourier transform of Ci = (cf,cf)^ 
and Dg = -eiEaa'Gaa'('?)(/,"'^>F", with F'^ = p + 
cos0„CT^ +sin0„(T3, and (/^'f) = Efc(4+g-^"cfc). 
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FIG. 3. Energy per site in CE (solid line) and C (dot- 
ted line) states as a function of ei/t and the corresponding 
charge disproportionation with (a) G'^-*, (b) G^^\ (c) G'"^^ 
and (d) G, with /3 = 0.5. In G'^> K[/Ki = 0.5, and in G 
K2 = 2K3 = OAKi. 

The full Hamiltonian H = H^^£ + Hl^^f can be solved 
by iteration. In the next we will make a comparison of 
the magnetic CE and C states. These two states have 
the same magnetic energy —Jaf per site, where Jaf is 
the AF magnetic superexchange between neighboring lo- 
cal spins, so that the relative stability of them is inde- 
pendent of the parameter Jaf- While other magnetic 
states, say, the FM and layered (A)-type AF states, have 
different magnetic energies so that their stabilities de- 
pend on JaFi and become less stable with respect to the 
CE and C state when Jaf increases. In the C state there 
is only one effective orbital in each site so that U has no 
effect and Hde — "^Sfe^^^c^ in the x-orientated FM 
chain. Such a property make the competition between 
the C and CE states alone very interesting. In the ab- 
sence of the electron-lattice interaction, without U the 
energy per site is E^^ = -0.693i and E^ = -0.637i, 
when U increases E'-^ keeps unchanged while E'-^^ ia.- 
creases and becomes higher than E'^ at about U = 5nl3, 
indicating that the strong electronic correlation would 
destabilize the CE phase towards the C state. When the 
electron-lattice interaction is taken into account. Fig. ^ 
shows the energy per site and the charge disproportion- 
ation dc — {n2j — n2j+i} as a function of ei/t in the 
CE and C states. Fig. y (a-d) corresponds to the ten- 
sor G^^\G^'^\ G^^^ and the present G, respectively. At 
e; = the energy E'-^^ > E'-' , the CE state is unsta- 
ble. When the electron-lattice interaction increases, it 
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is found that in the NLD case of G^^^ the CE phase 
always has higher energy than C, while in all the CLD 
cases of G'"'^', G'^' and G, there is a crossover from C 
to CE state with the increasing ci. So here NLD is not 
enough to stabilize the observed CE state, to obtain the 
CE state the cooperative nature of the LD must be taken 
into account. The different results in the NLD and CLD 
cases can be understood from Eq.(^. In the CLD cases, 
when a pair of neighboring sites are both occupied, the 
additional gi coupling favors different orbitals on the two 
sites. Since the orbitals on neighboring sites are the same 
in the C state but are different in the CE state, the lat- 
ter is more favored. At large charge dispropornation, one 
of the neighboring sites is almost empty and the gi cou- 
pling makes no difference between C and CE states, so 
in Fig. H (a) and (c) the energy difference 6E of the 
two states no longer increases with further increasing e;. 
On the other hand, 6E keeps increasing in (d), which is 
related to the g2 coupling in Eq.(Q), and its effect will 
be discussed later. In the calculation with G, it is also 
found that with different values of K2 and > we 
get qualitatively the same results. The calculated rela- 
tive displacements of the Z, O and Mn sites |Aq|/|<5q| 
and |^q|/|Mq.xl is independent of ei, the former actually 
depends only on K3/K2, and the latter decreases with 
increasing K2/K1 or decreasing K3/K1. In Ref. || the 
ratios | Aq|/|^q| « 0.56 and l^ql/l^q.^rl^ 0.93 were mea- 
sured at q = (7r/2, 7r/2, 0). In the present calculation 
with K2 = 2K3 = OAKi the two ratios are 0.57 and 
0.90, quite close to that in Ref. ||. 

At strong electron-phonon interaction (e/ ^ t), the 
charge disproportionation tends to 1, the electrons be- 
come localized and can be treated as classical objects. 
A classical treatment of electrons can simplify the study 
considerably and show clearly the effect of CLD, and in 
fact should be appropriate in some manganites in which 
the charge difference between neighboring Mn sites is 
close to 10. In the classical case, riq = ((5q,o+^q.Q|| )Vn /2 
in both the C and CE states, where Qy = (tt, 7r,0) 
and N is the total number of Mn sites, and mq — riq, 
= rriq = — nq/2 in the C state, mq'^ — riq/ A ± 

3(<5q,Q||/2 + (5q,_Q||/2)ViV/8,m^ = -nq/2 in the CE 
state. For GW(i = 1,2,3), the energies obtained from 
/q = Priq + rriq and Eq.(||) are the same in C and CE 
states. If we further take the Wigner crystal (WC) state 
into account, which has the same COSO as that of CE 
in the xy plane, but along the z direction the charge 
density is altering instead of stacking, and /q'^^ = 

[/«.CE(i^g.,.) + ^a^CE(l_g,,.)]/2 with = (0,0,^), 
then the energy difference between CE (or C) and WC 
states are ei{P-l/2)^K[/{2Ki+4:K[)- Jaf, -Jaf and 
ei{/3 - 1/2)^ / A~ Jaf for G^^), G^^) and G^^). So that in 
the cases of G^*) {i = 1,2,3), without Jaf WC should 
be more stable than CE state, and for /3 = 1/2 a finite 
Jaf would yield stable degenerate CE- and C- stacking 
states. On the other hand. Fig. ^(a) shows the energy of 
the C, CE, WC and stripe phase (SP) states with G at 



P = 1/2 and Jaf = 0, in which CE state has the lowest 
energy. A Monte Carlo (MC) simulation on 8 x 8 x 8 
lattice with periodic boundary conditions, in which we 
consider three possible electronic states on a Mn site in- 
cluding occupied by elongated orbitals d^^2_^2, d^y2_j,2 
and unoccupied, is performed in real space to find the 
charge and orbital configuration with the lowest energy. 
The compressed orbitals such as (1^2 _y2 ace- pot taken 
into account due to the anharmonic effectsli^'tj. For a 
given configuration, we calculate its /q in momentum 
space and get the energy through Eq.(^. Fig. ^(b) 
is the calculated phase diagram. Note that phase dia- 
gram is obtained not by comparing the several states in 
Fig. ^(a), instead, each state here has the lowest energy 
among all the possible charge and orbital configurations 
within the range of consideration. For /3 close to 1/2, the 
obtained COO is the same as that under CE-type spin 
environment shown in Fig. |l|(a). The striking feature of 
charge stacking (CS) along the z direction is pproduced. 
Such a stacking is usually attributed to Jaf", yet here 
our calculated results provides another possible explana- 
tion, that the CLD may also lead to the CS. It is worth 
mentioning that the above simulation can be easily gen- 
eralized to 2/3 doping, where the obtained COO for (3 
close to 1/2 is a CS state same as that observed in Ref. 
^, and such a CS can not be explained by Jaf- The COO 
at half-doping may be interpreted by the effective inter- 
action shown in Eq.(^. The COO in the xy plane can be 
explained by the gi and 32 terms. The gi term favors the 
charge ordering peaked at (7r,7r). In such an ordering, if 
site i is occupied, then i + a {a = x oi y) is empty and 
i -I- 2(5 is occupied. The 172 coupling between the occupied 
sites i and i + 2a then prefers the orbitals in the two sites 
to be different, thus the desired in-plane COO is formed. 
The stacking in the z-direction is related to the g'^ cou- 
pling which favors neighboring sites occupied by the same 
orbitals. When (3 is close to 1/2, along the z direction 
the gi and 32 coupling is effectively very weak as at an 
occupied site ff = (3 — 1/2 is small, then the g'l term is 
dominant. Note that here the COO is obtained without 
invoking magnetic interactions, so that the COO transi- 
tion temperature Tco can be higher than the magnetic 
transition temperature T/v , in agreement with the exper- 
iments in some doped manganitead Since the effect 
discussed above exists beyond the classical limit, in more 
general cases the CLD should also favor such a COO. 
Once this COO is built, the CE-type zigzag magnetic or- 
dering below T/v can be understood from the competition 
between the anisotropic electronic hopping and Jaf ■ For 
an occupied site i of orbital d;jj.2_j.2 [d3y2_r2), the elec- 
tronic hopping between sites i and i + x {i + y) leads 
the spins of these two sites to be parallel. On the other 
hand, the spins of i and its neighbors in the other two 
directions are antiparallel as the electronic hopping in- 
tegrals in these two directions are much smaller and not 
enough to overcome Jaf- In this way naturally the CE- 
type zigzag magnetic ordering shown in Fig. |l| (a) is 
obtained. In this picture of the COSO, appropriate at 
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least for those whose Tco > Tn, COO has its origin of 
the cooperative nature of the lattice distortion, and the 
spin ordering is the consequence of such a COO. 
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FIG. 4. In the classical treatment of electrons and with 
fixed Kz/K2 = 0.5, (a) engery per site in the CE, WC, C 
and SP states at /3 = 0.5 as a function of KijKx^ (b) phase 
diagram from MC simulation. The 4x4x2 unit cell of SP 
is shown in (b), where the two 4x4 lattices are in successive 
x-y planes, and x, y, o represent orbitals do,xi_^i, d^y2_^2 and 
a hole. 
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